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ABSTRACT 


A study has heen performed on accuracy aspects of stereo side-looking 
radar. The geometry of the radar stereo model and factors affecting visual radar 
stereo perception are reviewed. This is then employed to define limits to the 
vertical exaggeration factor of stereo radar and to compare it to that of camera 
photography. 

Radar stereo model accuracies 8tre analyzed with respect to coordinate errors 
caused by errors of radar sensor position and. of range, and with respect to errors 
of coordinate differences, i.e,, cross-track distances and height differences. 
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SECTION 1 


INTRODUCTION 


Stereo viewing of side-looking radar images can enhance the interpretation 
of the images by Improved observation of morphological details (Koopmans, 1973), 
it can be used to measure slopes and relative height differences (Dalke and McCoy, 
1968 ), and it can improve the accuracies of cartographic mapping and point posi- 
tioning (Gracie, et al., 1970; DBA-Systems, 197^+; Leberl, 1978). 

Generally the term stereo refers to visual perception when an observer 
views an overlapping stereo image pair and in his brain forms a three-dimensional 
replica of the imaged area. Stereo may also be used for a computation process 
employing monocular measurements in overlapping images: this may occur with data 

that generally do not permit visual stereo but may show horaologue details in 
overlapping images that can be identified monocular ly. 

Radar stereo is perceived like its photographic equivalent, although geom- 
etries are quite different (Figure 1). The advent of aircraft side-looking radar 
(SLR), more recently of satellite radar, e.g. , in the Apollo-17 and Seasat mis- 
sions, and several proposed satellite radar projects such as the Space Shuttle 
SIR-A, SIR-B, Venus Orbital Imaging Radar (VOIR), have led to a requirement to 
understand stereo radar. 

This report presents some basic considerations on stereo radar, addressing 
its visual perception and a geometric analysis of its possibilities and limita- 
tions. Essentially those aspects will be treated that are valid with synthetic 
aperture radar, since this is being used on spacecraft. 


HEIGHT 




Figixre 1. Comparison of Basic Geometries for Orthogonal, 
Central and Range Projection 


1-1 


SECTION 2 

VISUAL RADAR STEREO 


2.1 BINOCULAR VISION AND STEREOSCOPIC VIEWING 

Visual stereo perception has been analyzed by Aschenbrenner (1952), Fichter 
(195^*), LaPrade (1972, 1973), LaPrade, et al. (1975) and others. As a result, 
stereoscopic viewing is well understood, althou^ there remain unresolved 
questions . 

In binocular vision, the eyes center on a point P (Figure 2), so that a 
convergence angle y is subtended. A second point A generates images a"*, a'“* on 
the retinas, producing a so-called retinal disparity a that lesuis to observe the 
point A at a distance different from point P. According to LaPrade, et al. 
(1975), the retinal disparity angle must at least be 3" to be observable, and 
shovild be smaller than 1° for stereoscopic fusion of the two images. In Nowicki 
(1966) the minimum retinal disparity is given as 10" to 20" . 

Binocular vision can be simulated by presenting to the eyes two overlapping 
images of an object. We then speak of etereoeoopio viewing. We will not dis- 
cuss the various means or aids to stereoscopic viewing of overlapping images, 
nor the subjective problems that may arise in visueQ.ly estimating heights and 
slopes. These problems have been elaborately treated for camera photography by 
LaPrade (1973) and others. We will concentrate instead on stereo computations 
that follow the stereoscopic measurement of radar image coordinates. 
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Figm-e 2. Binocular Vision and Retinal Disparity a 
(From LaPrade et al, 1975) 
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2.2 RADAR STEREO ARRANGB94ENTS 


The most common radar stereo arrangements are illustrated in Figures 3a, b 
denoted by eaone^side and oppoeite-eide , and obtained in two separate flights 
or orbits. They were for the first time proposed by LaPreide (1963) and later 
analyzed ly Innes (196U), Rosenfield (1968), LaPrade (1970), Grade, et sd. 
(1970), Konecny (1972), DBA Systems (197*^), Goodyear (197U), Derenyi (1975) » 
Graham (1975» 1976), Leberl (1972, 1975, 1978). An analysis of visual stereo 
radar perception was attempted by LaPrade (1970), Graham (1975, 1976) and 
Leberl (1978). 

Apart from common same-side and opposite-side arrangements, stereo radar 
has also been proposed using images taken cross-vise (Leberl, 1972; Graham, 

1975), at different altitvides (Leberl, 1972), or in a single flight line or orbit 
using some sort of convergent scheme (Leberl, 1972; Carlson, 1973; Bair and 
Carlson, 197*^, 1975). These arrangements are illustrated in Figure U. 

Single flight line stereo cannot be realized with synthetic aperture rsdar 
(SAR). One may be tempted to believe that in a squinted mode SAR operation 
(non-zero Doppler correlation) one obtains two stereo looks by once imaging with 
a positive Doppler frequency (looking forward) and once with a negative frequency 
(looking backwards). It was pointed out that such arrangements would not pro- 
duce stereo (Leberl, 1972): Relief displacements would always be at an angle 

of 90° towards the nadir line, irrespective of the amount of Doppler frequency 
used for image correlation. As a result there will be zero parallax in squinted 
SAR image pairs and no valid stereo. An additional outline on this matter is 
presented in Appendix A. 


ANTENNA 1 ANTENNA 2 ANTENNA 1 ANTENNA 2 



Figure 3. Basic Stereo Radar Arrangements 
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The single-flight stereo schemes of Figiire U could be realized with real 
aperture radar, tilting the radar antennas around horizontal and vertical axes 
and using conical beams. Since these schemes are invalid with SAR they will not 
be pursued here. 


2.3 FACTORS AFFECTING VISUAL STEREO RADAR 

In order to visually perceive a three-dimensional model of the object, 
the two images of a stereo pair must be sufficiently similar: the image quality, 
object illumination, tones and textures must be comparable and the geometric 
differences (stereo parallaxes) must not exceed a value of 1° (retinal disparity). 
In aerial photography this rarely presents a problem, since sun angles do not 
change drastically in overlapping photos. However, illiunination in the active 
radar system depends on sensor position. Considerable differences of illumina- 
tion must therefore be anticipated in radar stereo pairs. 

Figures 5 to 8 present examples of radar stereo pairs demonstrating some 
of the 1; nits to stereo viewing. Figure 5 shows part of the Estrella mountain in 
Arizona (USA), imaged with an opposite side arrangement from an aircraft at 12 km 
altitude. It can be observed that slopes reflecting strongly in one image are in 
the radar shadow in the other image. A stereo impression can be achieved in the 
flat areas of this stereo pair but becomes impossible in the mountain. Figure 6 
demonstrates with a same-side stereo pair that viewing does not present any 
problems. Koopmans (1973) could show in some tropical areas of Colombia, that 
this type of same-side radar stereo produced a drainage analysis that was more 
complete than that from available (poor-quality) aerial photography. 

However, there exist also limits to successful same-side stereo. Fig- 
ures 7 and 8 present two Apollo 17 lunar satellite stereo pairs with same-side 
gewnetry and a very small stereo base. Since look angles are steeper than those 
of Figures 5 and 6, one obtains larger relief displacements and variation of 
image appearance, even with a small stereo base. In the flat parts of Fig- 
ure 7, stereo viewing is easy. However, in the Apennin mountains stereo fusion 
becomes nearly impossible due to the differences of image contents. In the 
rugged Oriental region of Figure 8 this is even more distinct. 

One may thus conclude that radar stereo viewing depends on: 

The stereo arrangement; 

The look angles off nadir; 

The stereo intersection angles; 

The noggedness of the terrain. 

In flat or gently rolling areas stereo will hardly present a problem. In 
rugged terrain, however, stereo is only possible with same-side geometry and 
improves with shallower look angles. For good visual stereo perception one would 
prefer small stereo intersection angles: this results in image pairs with little 

differences in tone and texture, however also in gecxnetry. For good topographic 
expression (vertical exaggeration) one requires large Intersection angles. There 
exists thus a trade-off between geometric accuracy and ease of perception. 
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Aircraft Stereo Radar Ima^-'e Pair, Saii.e Side Gefimetry 
X-Band, IP-lun Altitude; Estrella Mountains, Arizona 
(Courtesy of Aeroservice Goodyear) 


Figure 6 
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For the stereo observation of man-made objects, LaPrade (1975) reports 
optimum stereo viewing if look angles vary between 37® and 67® off nadir, with 
intersection angles of about 12® to 15®. 
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SECTION 3 


GE(»4ETRY OF THE RADAR STEREO MODEL 


3.1 GROUND RANGE PRESENTATION 

Radar projection lines (lines connecting image and object) are circles that 
are concentric with respect to the flight line (Figvire l). 

If all observed slant ranges are projected on a datvan plane then the radar 
image is in gvound range presentation (Figure 9). Relief displacement in this 
case is shown in Figure 10. 

Stereo imaging depends on parallaxee. One defines abeotute parallax as 
shown in Figure 11: the two radar images of a stereo pair are superimposed 

onto each other in such a way that the flight lines coincide. Absolute parallax 
then is the distance between the two homologue images A', A." of an object point 
A, If an area is flat, and if there is ground range presentation, then all 
points have the same absolute parallax — the stereo model appears flat. 


Figxires 12 and 13 show differences of absolute parallax. These parallax 
differences may, for example, be measured using a parallax bar and stereoscope, 
and must then be converted to terrain height differences h. Expressions for this 
were formulated by Koopnans (197^) » Derenyi (1975) » Leberl (1975) and others. 
These expressions use often an approximation by replacing projection circles 
(wave fronts) by straight lines normal to the line of sight connecting object and 
antenna. The effect of neglecting the curvature of the projection circle reduces 
with the distance from the antenna and is thus very small for satellite radar. 


Formulas to compute height h from radar parallaxes Ap are then, according 
to Figures 12, 13; 


Same Side 


Opposite Side 


p'' = h • cot 6" 


p" = h • cot. O' 


p" = h • cot e" • 

p" - p' = Ap = h (cot 6" - cot 0') 

h » Ap/(cot d" - cot 0') 


p"' = h • cot 6" 

p' + p" = Ap=h (cot 0" + cot 0') 
h = Ap/(cot 0" + cot 0") 


( 1 ) 


Practical application of formulas (l) depends on the knowledge of look 
angles 0', 0'*''. We see that these vary for each image point. Also, they are 
only approximately known, since not the angle to point P is available, but the 
angles to points P"*, P". 
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Figure 9. Radar Image Presentations 
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Figure 13. Stereo Radar Parallaxes, Opposite Side Geometry. Shown is 
also the camera arrangement to obtain the same parallax 
(LaPrade, 1970). 


Expressions (l) are the radar equivalent of the simple photogrammetric 
parallax formula: 


h = Ap • H/B 


as used in photointerpretation. 

It may be significant to point out that in photography object heights h 
will be directly proportional to parallax difference Ap, involving merely a con- 
stant proportionality factor H/B. With radar, however, parallax differences Ap 
produce different heights h, depending on the look angles 9', Q". 


3.2 SLANT RANGE PRESENTATION 

A single slant range radar presentation generates a panoramic visual impres 
sion as illustrated in Figure II 4 and confirmed in the images of Figures 8 and 9. 
Also a stereo-image pair creates a bowed stereo model. For slant range 



Figure lU. Visual Impression Obtained 


APPAHENT-^ 

TERRAIN 

SURFACE 


Slant Range Presentation 


presentations, one seems to have to deal with a more complex situation, since 
the absolute parallax is veuriable even for points in the datum plane. 

The absolute parallax p must be defined as: 


P 


r - r 


(2) 


where r", r" are slant ranges to etn object point, measured in images (l), (2), 
respectively. From Figure 15 we know that: 


- ^ 2 . 1/2 

r = (H y ) 


r" = (H^ + (y-B)^)^/^ 


so that 


P 


(H^ + - (H^ + 




(3) 


The absolute parallax depends on the value of y. It can be shown that 
the relationship between p and y defines a super-hyperbolia cylinder in the 
x»y»P “ space (Leberl, 1978). One can verify, that for y = B/2, p = 0; for 
y > B/2, p > 0; for y < B/2, p < 0. These facts are illustrated in Figure 15 
and clearly present the super-hyperbolic cylinder of flat terrain. Relief is on 
top of the cylinder. According to Figure l6 the height h above the super- 
hyperbolic cylinder results from: 


p' = h • cos 0' 
p" = h • cos 0" 

h = Ap/(cos 0' + cos 0"’) (U) 

Where the minus sign is for same-side, the plus sign for opposite-side 
arrangements . 


3.3 RIGOROUS STEREO INTERSECTION 

Rigorous radar stereo intersection algorithms were used by Gracie, et al. 
(1970), DBA-Systems (197^^), Leberl (1976), and for underwater SONAR images by 
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Figure 15. Stereo in Slant Range Presentation, Terrain Rolled 
onto Cylinder 


O' 0" 
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Clerici suid Konecny (1978). These formulations define a position vector £ for 
an object point P, using all available information on sensor position £, attitude 
s_ (velocity vector), squint angle t, slant range r. Figure 17 illustrates the 
situation: 



l£ - ^"\ 

= r" 



(b) 

(£ - s') 

- sin t( ' 

\i'\ 

l£ - s.' I) 

(c) 

(£ - S") 

= sin t( 

ir- 

! l£ - s"| ) 

(d) 



Figure 17. Rigorous Stereo Intersection 
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where Equation (6) presents the two range spheres (a) and (b) and Doppler 
cones (c) and (d). The elements of vector £ are imknown. Since U equa- 
tions exist for 3 unknowns, an adjustment can be carried out. The algorithm 
for this is listed in Appendix B. 
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SECTION k 

VERTICAL EXAGGERATION IN STEREO RADAR 


LaPrade (1970) and Graham (1976) have discussed vertical exaggeration in 
radar stereo models. These considerations relate the radar stereo parallax to 
the base to height ratio Bc/H of equivalent camera photography as shown in Fig- 
ures 12, 13. From these figures it is obvious that for ground range presenta- 
tions (compare Equation (U) with B^/H ■ dp/h): 

B /H » cot 9' - cot 9" (7) 

c 

and for slemt range presentation (compare Equation (U)}; 


B /H = cos 9' - cos 9" (8) 

c 


where B^ is the stereo base of the equivalent camera arrangement. 

The ratio Bc/H is a measure of vertical exaggeration q that has been 
extensively analyzed in photogrammetry. The factor q can be defined as the 
ratio between the true height-to-base-ratio h/b of an object and its visually 
perceived stereo model equivalent, 

q « (h • bj^)/(hj^ • b) (9) 


Following LeiPrade (1972), the factor q for photography is found to be; 


q 5 • B /H (10) 

c 

if the stereo-model is observed in a stereoscope and if the separation of 
photographs produces an angle V of 0.20 radians (y as defined in Figure 12). 

One obtains then as the vertical exaggeration of stereo radar: 


q a 5 (cot 0' - cot 6") (11 ) 

or, as shown by Graham (1976) and deriving from Figure 12; 

q « 5 * H • B/{y(y - B)} (12) 
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Figure 18 presents the values for q for airborne radar and the proposed VOIR 
orbital mission. For the optimum stereo radai* arrangement as defined ^ 
LaPrade (1975)» we find (compare Section ?.?); 


37 ® < S'* < 67® 

12® < A0 < 15® 


where 


A 0 = 6' - e" . 


Therefore 


2.3 > q > 1.8 . 


These values for q are smaller than what one would obtain with wide-angle 
photography, where B/H = 0.6, so that: 


q = 3 . 


It is also obvious that q is variable across a radar stereo model while for 
photography it is a constant for a given stereo pair. 

For VOIR, vertical exaggeration must be expected to only amount to 
q 0.6, since stereo intersection angles are small: 3?3 > A0 > 2?8. 


3.3® 3.0® 2.rf* 



365 395 425 y(KM) 


ir q for Aircraft and VOIR-Stereo Radar, 
near to far range 


SECTION 5 


ERRORS OF THE RADAR STEREO MODEL 


Errors of the radar stereo model shoiild he divided into two groups: 
absolute errors are errors of individual point coordinates, or planetocentric, 
positions and height. Relative errors are those of coordinate differences.* 

Error analyses have been previously published by Rosenfield (1968), LaPrade 
(1970) and Leberl (1972) for airborne radar, and by Leberl (1978) for satellite 
radar. These analyses were limited, e.g., to compare same-side with opposite- 
side stereo, and addressed only absolute errors. The present analysis aims at 
a more general approach, identifying both absolute and relative errors irrespec- 
tive of flight arrangement. 

Figure 19a defines absolute errors that result from errors of the orbit 
positions s^'i of the velocity vectors of slant ranges r", r" , 

while Figure 19b shows how a pair of points with absolute errors combines to 
define relative errors. 

One is dealing with error vectors , ds^. ds^^ and with dr", dr"'. 


5.1 COORDINATE ERRORS 

For the derivation and discussion of errors of the radar stereo model, a 
coordinate system and flight arrangement as shown in Figure 20 are used, where: 


» (x^', y^', z^') = (0, 0, H) 

i' = (i^', 0, 0) 

s" = (x^", y^". z^") = (0, B, A) 

s" = (x^", 0, 0) 


(13) 


•Relative errors are often used as defined here. Correctly, however, this term 
should be used only for the ratio of error and observed quantity and be 
expressed in percent (%) . 




Figure 19. (a) Absolute Errors of Stereo Model 





FUGHT LINE 



Figure 20. Coordinate System 


It is obvious to see that an error dx^" has no effect on y and z, so that; 


dx " 

= dx '12, 

dx ^ 

= dx 

'72, 


xo 

o 

xo 

o 


dy 

= 0 

dy 

•'xo 

" = 0 


(iM 

dz ' 

= 0 , 

dz 

" = 0 

• 


xo 


xo 




The effects of dy^,', dyQ^ 

are less 

obvious . 

Figure 21 

clarifies 

the relation- 


ship. There is no effect on x. In y, we find: 


dy H 

I ~ r' 
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Figure 21. Errors Due to -Component of Sensor Position 


with 


JL - sin e" 

sin (O' - d") 

So that : 

dy = H • sin 6"/(r' sin (O' - 0'^)) (15) 


But: 



sin (90° - 9") 
' sin (0' - 0") 


r' sin (0' - 0") 


B sin (90° - 0") 


(16) 


Introducing (l6) in (15) gives: 


dy = dy^" H tan 6'7B 


so that, finally; 


dyy^" = (B-y)/B (17) 


We can substitute for y/B also; 


z _ 

y 

. L' 

sin 

6' sin 

(90° + 0") 1 

B 

r 

’ B 


sin (6' 

- 0") ■ 1 - cot d' tan e" 

A similar development 

can be 

found for dz ", so that we obtain: 
yo 

dx 

yo 

= 

0 



• 

dy 

yo 

= 


" (B - 

y)/B = 

dy + cot 0" tan e") 

0 

dz 

yo 

s 


' y(y - 

■ B)/H B 

= -dy^"’/(cot 6" - cot 0') 

and similarly for the 

other sensor station: 

dx ' 
yo 

= 

0 




dy 

yo 

= 

dyo' 

• y/B 

= dy^ 

"/(I - cot £J" tan U") 

dz ' 

yo 

= 


' • y • 

(y - b)/(h • B) = dy^V(cot e" - cot e') 


(18) 


Equations (l8) 
linearly; this 
representing a 
.o the y^-term 


show that the y-diraension of the stereo model is deformed 
is a scale error. The model height is deformed linearly in y, 
model tilt, and in addition there is a nonlinear model bow due 
Error curves are presented in Appendix C. 


The effects of height' errors dZg' , dzQ^" are shown in Figure 22. A 
derivation like the one for Equation l8 results in 
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A 

i 


< 



dx ' 
zo 

• 

0 

dy 

zo 

= 

-dz " H/B = dz '/(tan 0" - tan 0') 
o o 

dz 

zo 


-dz ' (y - B)/B = dz '/(l - tan 0' cot 0") 
o o 

dx 

zo 

= 

0 

N 

O 

= 

dz " H/B » dz '/(tan 0' - tan 0") 
0 o 

dz '' 
zo 

= 

dz " Y/B * dz '/(I - cot 0' tan 0") 

O 0 


( 19 ) 


An error of sensor altitude does not deform the model y-dimension, but 
merely shifts the model in y. In height, there is a linear effect of model tilt. 
Error curves are shown in Appendix C. 
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Erroneous sensor attitudes, caused lay errors of the velocity vectors 
lead to nonintersecting projection circles. This is illustrated in Fig- 
ure 23. for the example of an error d^Q' which creates aji antenna swing k. We 
need to define a model point, that most logically is chosen halfway between the 
error- free point obtained if no dyo" error had been present. 

From Figure 23 we find that an error dy^' affects the x-coordinates of the 
model point P. shifting it to position p. but leaves the y.z-coordinates 
unchanged, at least to a first-order approximation: 


* dy^" • y/2 = dy^' • H • tan 0 72, 


dy. " = 0, 

•'yo * 


dz. ' = 0, 

yo 


® ^o^^ ^ ®"/2. 


dy. = 0. 

yo 


( 20 ) 


dz. = 0. 

yo 


TOP VIEW 



Figure 23. Error Due to dy^ -Component of Velocity Vector 
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Similarly, the dz ■*, dz '"’-components cause a tilt i of the radar antenna, 
so that: ° ° 


dx. ' 

zo 


dy. ' 

zo 


dz. ' 
zo 


dx. " 
zo 


<3y 


zo 


dz. 

zo 


dz ' • H/2, 

o * 


0 . 

0 . 



0 , 

0 . 


( 21 ) 


An erroneous dx^', dxQ "-component will not affect the antenna attitude and 
thus not produce errors of the stereo model. It will, of course, affect antenna 
position when integrated over time. This error was considered in Equation (lU). 

A last error of the stereo-model results from erroneous range data r', r". 
Figure 2U shows that: 


dx ' = 0 

r 



dr' • r'/B 


dr' (y^ + 


dr' • cos e'Vsin (6' - 6") 


dz^' = dr' • r' • (y - B)/B • H = dr' sin 0"/sin (0' - 0") 

dx " = 0 

r 


(22) 


dy^" = -dr" • r"/B - -dr" • cos 0'/sin (0' - 0") 


dz^" = -dr" • y • r"/BH = -dr" sin 0'/sin (0' - 0") 


These errors are nonlinear in the y and z-dimensions . Error curves are 
plotted in Appendix C. 
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Figure 2^. Error Due to Range Error 


5.2 ERRORS OF COORDINATE DIFFERENCES 

Coordinate errors generally are of lesser importance than are errors of 
coordinate differences — these are actually the entities that are to be derived 
from stereo radar. In this context it is legitimate to study only the effects of 
an error of dB of the stereo base B, of erroneous height differences dH between 
the two antenna posit i(5ns, and of range errors dr', dr". The along track model 
dimension is not affected by dB, dH, dr', dr" (compare Equation l8, 19, 20), 
Therefore, it is of little concern in a stereo-radar analysis. 


5.2,1 Stereo Base ; dB 

The errors of coordinate differences in cross track and height are denoted 
by dAy, dAz. They follow from partially differentiating the following equations 
for the coordinate errors ^ using dB = dy^" - dy^' (compare Eqn. (l8)): 


dyg = dB + const. 


( 23 ) 
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The partial derivative with respect to y and z produces the following 
equation for the errors dAy, dAz of coordinate differences; 


(2l*a) 


(24b) 


Elquations (24) describe the change of coordinate errors with varying the model 
coordinates y and z. It is to be noted that the model height z and flight height H 
complement each other, so that with the definitions of Figure 20 dH = dz and 
AH » Az. 


Equations (24) are valid for small values of Ay, Az — they refer to errors 
defoming small coordinate differences Ay, Az in model space. Conclusions to be 
drawn from Equations (24) are evident; 

• Cross-track distances have errors directly proportional to 
base-errors . 

• Height differences have a more complex error behaviour; for vertical 
structures (Ay = 0) errors increase with both a linear and a non- 
linear term in y. 

• Heights measured of points that are spaced apart for Ay have errors 
that consist of a small constant terra and a large term growing 
linearly with y. 

Numerical examples are presented in Appendix C. 


Ay 

- -r i' 


I Ay - ~ Az ! dB 
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These equations, vhen partially differentiated, lead to: 


dAyjj = — dz ( 26 a) 

dAzR = “ dy ( 26 b) 

We find here a slxtq)le error situation: cross track distances and height dif- 

ferences have errors that are directly proportional to dH/B, trhere cross track 
distances have errors only if end points are at different altitudes, and heights 
of vertical structures would be error- free. 


5 . 2.3 Range Errors: dr*, dr"^ 

Range errors may be of two origins: 

(a) Limited range resolution causes random errors of range measurement. 

(b) Errors of range calibration, measurement etc. may lead to systematic 
range deformations. 

Range resolution defines the height definition ^ thus the acciuracy with 
which we can set a mark on the surface of the observed stereo model. Similarly, 
also the y-coordinate is subject to a limited definition. 

An expression is easily found for the cross-track and height definition 
starting from Equations (22), where the random errors of height and cross-track 
coordinate result from error propagation and are denoted by: 


cos^ Q' + cos^ e''P^/sin (0' - 9") 

(27a) 

sin^ 0' + sin^ 0"|^^^/sin (0' - 0") 

( 27 b) 


or: 

Oy = Oj. jr'^ + ( 27 c) 

0=0 |r'^ (y - B)^ + • y^P'^^/BH (27d) 

z r ( ) 
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An actual stereo-model warp may result from systematic range errors. Cross-track 
distances and height differences would be affected as follows: 


ddy^ 


dAz 

r 


/dr' y dr" (y-B) \ . . /dr' • H dr" • H\. v 

(r7^ - -) - B r" P 


/dr' • r' , dr' (y-B) • y dr" • r" dr" • y • (y-B) \ ^ 

\ B • H * ~ B . i n r*" ■ E • It i . It • r” * 

♦ I (y~B) _ dr' * ^ _ dr" » y dr" * 


^ * y r \ 

B^ ? B • r" B • I? f 


(28b) 


From these algebraic expressions we cmi conclude that warping of the stereo model 
is nonlinear. However, all error coefficients are smaller them 1 so that model 
height differences and distances can be expected to have errors smaller than 
those of the source. Again, numerical values eu-e presented in Appendix C, 

5.3 DISCUSSION 

From f'e algebraic expressions derived in the previous sections it follows 
that different results are obtained for errors of coordinates versus distances or 
height-differences. The formulas and numerical evaluations in Appendix C 
(Figures Cl, C2) show that coordinate errors will be large in both aircraft and 
spewecraft radar stereo: for aircraft, with comparatively large stereo inter- 

section angles, coordinate errors will be 1 to 5 times larger than the error of 
sensor position or of range measurement. Heights z will in the near range be 
more accurate than the across-track coordinate y; at far ranges the situation is 
reversed; y is more accurate than z. 

For satellite radar such as VOIR — without two distinctly different look 
angles for overlapping swaths — the coordinate error will be 8 to 15 times 
larger than the error of sensor position or of slant range. 

In many instances, however, it is not the absolute accuracy that is 
important, but rather the error of a radargrammetrically measured distance 
and height difference. 

As illustrated in Tables Cl, C2 of Appendix C, both across-track distances 
Ay and height differences Az are less sensitive to errors of sensor position 
(stereo base B, sensor height H) or of range r than was found for the coordinates 
themselves . 

In comparing we find, that at the center of an aircraft stereo model such as 
that shown in Appendix C an error dB of the stereo base B of 100 ra will change 
the observed height of an object point for 26l m, its across-track position for 
237 m. Were one to employ VOIR-stereo, then the sensor error dB » 100 m would 


f 

J 

i 

result in a height change of 935 m, and a cross-track change of 988 m. A 
height difference, measured across a horizontal distance Ay would be affected 
differently in the same area of the stereo model: a base error dB * 100 m would I 

result in an error dAz of the height difference As = 1 km (Ay « 1 km) that ! 

amoiints to -11 m (aircraft) or -2.5 m (VOIR). I 

i 

The error propagations are similar when considering the effects of other 
error sources. Systematic errors of slant range have the smallest effect on 
height differences or distances. However, limited range resolution does not lead 
to systematic but to random errors. These generate an inaccuracy of measuring 
height differences and distances that are a factor larger then the error of a ! 

single coordinate meastu’ement . Thus a range resolution limited to +100 m for 
VOIR will create random errors of height differences of ±1.1 • ^ to^ 

±1.1* • /2 km, and random errors of cross track distances of ±1.3 • to 
±1.1* • km. 
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SECTION 6 
CONCLUSION 


We have reported basic considerations on stereo side-looking radar, 
presented a general algorithm for stereo radar intersection and derived algebredc 
expressions for the errors of point coordinates obtained by stereo-radargrammetry. 
We have, based on these formulas, given expressions for deformations of stereo 
models. These deformations concern cross-track distances and height differences 
and may thus be called "relative accuracies" of coordinate differences. The 
formulas are presented as derivations and essentleO-ly relate to systematic errors 
of the stereobase, B; of sensor height, H; and of slant ranges, r", r". 

The algebraic expressions vere derived using simple geometric considera- 
tions. They permit a wealth of conclusions to be drawn on stereo-radar. Some 
numerical evaluations of the formulas are presented in Appendix C for a general 
aircraft stereo radar arrangement and for a proposed satellite project VOIR. 

Radar stereo can be evaluated using vertical exaggeration as a figiire of 
merit. It was found that the vertical exaggeration is limited for satellite 
radar stereo if no distinctly different look angles are employed. For VOIR, the 
vertical exaggeration would not go beyond 0.6, while for aircraft radar optimum 
visual stereo would produce an exaggeration factor between 1.8 and 2.3. For 
camera photography, this same factor usually amoimts to 3. 

Errors of coordinates and of distances meas\a*ed in the radar stereo-model 
are systematic when caused by slowly varying errors of the stereo base or of 
the sensor height. Coordinate errors may be 10 times the size of errors of 
height differences when using aircraft radar; for VOIR this factor may amount to 
several hundred; this means that systematic errors of the stereo base, or of the 
sensor altitude, have a very slowly varying effect on the radar stereo model. 
Random errors, however, such as those resulting from limited raden: resolution, 
have an effect on distances and height differences that is essentially of the 
same magnitude as the effect on coordinate errors. 

Presently there is only limited knowledge on the optimum radar stereo 
arrangement; stereo intersection angles should be large from a geometric point of 
view, but small in order to permit visual stereo perception. Present indications 
are that visual radar stereo is possible with intersection angles of up to 15®, 
where angles off nadir should be not smaller than 1<0*. These are angles that 
should be used in present planning for future satellite radar missions. Limits 
to range resclution will considerably magnify the errors of observed height 
differences and cross-track distances if stereo-intersection angles are chosen 
that are small. However, in order to understand the exact limits of visual stereo 
it is necessary that extensive experimental work still be carried out to analyze 
the Interaction of stereo-intersection angles, look angles off nadir, and terrain 
ruggedness in their effect on the subjective experience of visually measuring and 
interpreting a radar stereo model. 
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APPENDIX A 


THE IMPOSSIBILITY OP SINGLE FLIGHT STEREO SAR USING 
NON-ZERO DOPPLER PROCESSING 


It will be shown that processing of a single SAR image atrip with 
differential Doppler frequencies does not lead to valid stereo. 


Al. THE DOPPLER CONE AND RANGE SPHERE 

Imaging with SAR using a non-zero Doppler frequency for correlation is 
equivalent to imeiging with a (synthetic) radar beam of conical shape: from a 

given sensor position points will be imaged that produce the same Doppler fre- 
quency shift ~ this implies that they are on the surface of a cone (ccanpjare 
Figure A-l). The locus of an image point for which a slant range r from the sensor 
has been measured is a sphere, the so-called range sphere (compare Figure A-2). 


A2. DESCRIPTIVE GEOMETRIC CONSIDERATIONS CONCERNING RADAR PROJECTION LINES 
AND STEREO 

Given is a point P in a radar image and the sensor's position and attitude 
when the point was imaged. From the single radar record the geometric locus of 
the point P can be established as an intersection of the range sphere and 
Doppler cone. 



Figure A-l, The Concept of Imaging with "Squint” 
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It is trivial to see that the locus of point P is a circle concentric with 
respect to the flight line and in a plane perpendicular to it (compare 
Figure A-3). 

From this we can conclude that radar projection lines are circles. This 
fact was illustrated by Figure 9 in the main part of the report. 

If we assume a reference datum for the imaged surface, then we can define a 
" ground range preeentation'' according to Figure 9* Points outside this reference 
will show relief’-dieplaaement according to Figure 10. Radar stereo is based on 
different relief displacements in overlapping stereo image pairs, for example as 
shown in Figure 12 for stereo with two flight lines. It is obvious that no stereo 
is obtained with images that display identical relief displacements. 

From the above we can conclude that for a given single flight line one will 
obtain for any given object point exactly the same projection circle, whatever the 
Doppler cone angle t may be: two correlations with different Doppler frequencies 

(positive and negative, for example) will result in the same projection circle. 

An intersection is not defined (compare Figure A-U) . 


A-2 




Figure A-3. Definition of the Intersection of Range Sphere and Doppler 
Cone. Result is a "projection line" for an object point P. 


In other words: relief displacement will always be the same irrespective of 

Doppler frequency. Mountains will always fall towards the nadir line in a direc- 
tion perpendicular to that line. Therefore two images will not produce the desired 
stereo model of the object. 

However, one may obtain the visual Impression of an overall bow of the 
imaged area, if two differential correlations are viewed under a stereoscope. 

This may be due to an overall deformation of systematic nature, depending on how 
the correlations are being produced. That, however, is not a stereo model of the 
terrain relief. 
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Figure A-h. Intersection of Two Projection Lines in Single Flight 
"stereo" Arrangement with Non-2Lero Doppler Processing 
not Defined. Intersection Lines of Range Sphere and 
Doppler Cone Coincide. 


A3. ALGEBRAIC CONSIDERATIONS 

The facts of Section A2 will now be discussed algebraically. For clarity 
some assumptions are made about the coordinate systems. 

A coordinate system is defined with axes x» y, z, where x is parallel to 
the flight path. The sensor moves along a straight line in the xz coordinate 
plane (compare Figure A-5). 
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S(X .o,zj 

w o « 


NADIR LINE 


f(X,Y,2) 


Figure A-5. Definitions for the Formulation of Equations (Al) to (a6) 

We define the position of the sensor by the vector s* = (x'» 0, z ), the 
position of the object point P by the vector £ = (x, y, zT. The°cone angle is 
designated by t'. 

We find the equation of the range sphere as: 

(x - + (z - Z-)^ * U 

o o 

We find the equation of the Doppler cone as: 


(x - x') / sin j' » r' 
o 


The intersection of range sphere (Al) and cone (A2) is found by sub- 
stitution of (2) into (1); 


(A3) 


y + {z - z^) » r - r sin t 

o 

and 

* (z - cos^ t' 

This Is a circle in a plane x = constant. The center of the circle is at 
(x, 0, z^), its radius is r' cos x' 

Equation (A3) represents the projection line obtained from sensor posi- 
tion £" using a cone angle x'. 

A second correlation with a different cone angle x", but frcxn the same 
object £ = (x, y, z) corresponds to a sensor position * (x 0, z ), 
where x" results from Equation (A2): ° ° 

(x - x") / sin x" = r" (A^t) 

o 

The projection circle is given by: 

2 / \2 >>2 2 fnc\ 

y + (z - z) = r cos t (A5) 

o 

The second projection circle is also in the plane x * constant, with the center 
at (x, 0, z ). Since it has to pass through the point £ we obtain the 
identity: 


..2 


2 ^ 
cos T 


= r 




cos 


(a6) 


This identity confirms that both images produce the same projection circle, 
so that there is no intersection defined. We must conclude that no valid stereo 
is possible. 

Al*. CONCLUSIONS 

This appendix proves that single flight stereo with SAR is not possible on 
the basis of two correlations using different Doppler frequencies. 

A strictly rectilinear flight line was assumed. However, conclusions 
remain xmaltered with small deviations of the flights from a straight line. 
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AH>BHDIX B 

STRICT SOLUTION TO THE RADAR STEREO INTERSECTION PROBLEM 
B.l ALCK^ilTHII 

stereo radargraaoBetric computation may start out from Equations (6) 
of Section 3. If we assume that ve deal with SAR and that image presentation 
is without the effect of squint angle x (x = o), we obtain for each point: 


with 


£ - s'l 


|£ “ 


L' ' (£ “ iL"") “ 0 


i" • (£ - £") = 0 


2 = (x, y, z)' 


s. = (x^, y^, z ) 

o o o 


(B.l) 


i K’ 


This is a nonlinear system of equations for the unknown coordinates £ = 

(x» y, z)^. For a least squares algorithm we need to linearize: 

C.v + Dd£+f»0 (B.2) 


where d£ is the vector of corrections to the approximate solution £, v is the 
vector of observational errors, f is the vector of contradictions. We find; 


A£ = (Ax, Ay, Az)‘ 


V 't V V ", v. ", V. ", V. ", V ", V '" , V "" , 

xo yo zo xo yo zo r * xo yo 

T 

V "", V. V, "", V. V "") 

zo xo yo zo r 
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C? A-t 

wkedwhr,^ 


lb 

80 that an in^roved approximation £ of £ is found: 


£ » £ + d£ 

Matrix C is: 



S 

-2 (x - x^-) 

^1.2 

S 

-2 (y - 7 q' ) 

^1.3 

s 

-2 (z - 2^-) 

^1.7 

= 

-2 r' 

^2,8 

s 

-2 (x - X 

o 

^2.9 

= 

“2 (y - y^" 

^2,10 

= 

-2 (z - z 

o 

^2, lit 

= 

-2 r'" 

S,1 

= 

-X ' 
o 

^3,2 

s 


^3,3 

s: 

-z ' 

o 

S.l* 

= 

(x - x^-) 

^3,5 

X 

(y - y„') 

S.6 

s 

(z - z^') 

^U,8 

- 

-X - 

o 


8 


^4.10 

* 

-z '' 

o 

^4,11 

8 

(x-yo") 

"l:.12 

8 

(y - y^") 

‘^4.13 

8 

(z - z^") 


Hot FiUlED 

(B.3) 



All other elemnts of G are zero* For D: 



X 

2 (x - 

v> 

”1.2 

= 

2 (y - 

v> 

“1.3 

= 

2 (z - 


“2,1 

= 

2 (x - 

X ) 

0 

”2.2 

= 

2 (y - 

V'> 

®2,3 

= 

2 (z - 

V'> 

”3.1 

= 

X ' 

0 


°3,2 

= 

#o' 


”3,3 

= 

i - 

0 


'’-..I 

X 

X - 

0 


®1.,2 

= 

^o" 


®U.3 

X 

i " 

0 



For f: 


f 

f 


= - II- 8'j - r' 
fg = - |i- 

^3 * -L' ' 

= -s" * (i-s'O 


The least squares method leads then to; 

Ae = - (C g D)"^ (C g f 


(B.U) 
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where ^ is the variance-covariance matrix of the observed quantities x 


®o ♦ *o * ^o * * 

available a priori. 


• V V 

f • 


yo 


and mu§t be® 


The algorithm is iterative. After a preliminary computation to find an 
approximate solution £, one computes C, D, f_ and Ag.. With an improved 
^ s £ -f A£ one repeats the procedure until A£ becomes Insignificantly small. 


B.2 COMMENTS 

The above algorithm is simple and straightforward. In the event that not 
SAR, but resd aperture is available, one must substitute for £ the vector in 
the direction of the physical antenna. This direction can be described by 
pitch ^ and yaw angle k: 


X = cos 6 cos K 

o 

y^ = cos ^ sin tc (B.5) 

z « - sin ♦ 

o 

A totally different stereo computation results if Equation (5) of Section 3 is 
used as follows: 


s' + r' - s" - r" = 0 


(B.6) 


This is a coplanarity condition for the two orbit stations ('), (") and the 
unknown point. It consists of three equations with two unknowns fl", fl" since; 


= A r “ 
o 


fB.7) 


where 


r ' 
—o 


(sin T, (sin*^ ft' - sin^ 




COS 




(B.8) 


A rotation matrix 


Vector r^" gives the point in the antenna coordinate system, while A rotates 
that system into the reference coordinate system. 


B— U 


Details about the stereo computation with Equation (B.6) can be found 
In the literature (Leberl, 1976 a, 1976 b, 1978) together with a numerical 
comparison of the two alternatives. The approach with Equation (B.l) appeared 
to be somewhat superior In an application to lunar Apollo 17 - ALSE - VHP 
imagery with small stereo Intersection angles of 3*. 
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APPEUDIX C 

DEFOBMATIONS OF RADAR STEREO MODELS 


Equations (l6), (I9)« (22) » and (2M), (26), (28) have been evaluated for 
tvo stereo arx^ngements as showi in Figure C-1. For the aircraft radar, the 
ai^es of stereo intersection are 26** > Ad > 8?6, with 32?U < 8" < 68?2 and 
26?6 < < 59?5. For VOIR, 3?3 > A0 > 2?8, UU?2 < 6' < U8?6, U0?9 < 9" < ^5? 

Overlap of images is 60 %, 

Presented are error coefficients, not errors themselves. Figure C-2 shows 
aircraft stereo model coordinate errors dy, dz due to dyg'"’, &Zq" ^ dr', whereby 
the y-coordinate is plotted on the’ abscissa. For the example of dy^'*", a model 
error dy will be produced in the near range that is twice to 3 times the amount of 
dyg". It is obvious that the height error dz varies more than dy with increasing 
range, and that it is particularly affected at far ranges by dyo", dr'. 

Tables C-1, C-2 present the errors of height differences Az and distances 
Ay, again in the form of coefficients: for the exanq>le of a base error dB in 

VOIR, the height difference Az » 3 ion, measured across a Ay of 1 km at near 
range (y » 365 km) will have an error dAz of 0.017 • dB. 


• -40KM 




Figure C-1. Stereo Arrangements for Numerical Analysis 
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Figure C-3. Coefficients of Sensor Position Errors <Jyo'» and of Range Error dr' in the Equation for 

the Resulting Coordinate Errors dy, dz in the Stereo Model from VOIR, Figure C-1. 




Table C-1. Coefficients of Errors dB, dH, dr In the Equations for the 
Errors dAz of Height Differences Az, and for Errors dAy of 
Distances Ay, in an Aircraft Stereo Model of Figure C-1. 


Error Coefficients 


Base 

dB 

Height dH 

Range dr 

(km) 

Ay 

Az 

Ay 

Az 

Ay 

Az 

Az 

Ay 

y 



0.000 


0.000 

■1 

0 




0.081 

0.125 


-0.036 


1 

0 



0.2l*l* 

0.37; 


-0.107 

K8il 

■ 3 



-0.125 

-0.225 

IH 

0.125 

0.01*3 

0.01*2 

0 



-0.125 

-o.ii*i* 


0.125 

0.008 

o.oii* 

1 

1 

13 

-0.125 

0.019 

MM 

0.125 

-O.O6I* 

-0.01*2 

3 



-0.625 

-1.125 


0.625 

0.216 

0.211 

0 



-0.625 

-1.0l*U 


0.625 

0.180 

0.183 

1 

5 


-0.625 

-0.881 

m 

0.625 

0.109 

0.126 

3 



0.000 


0.000 

0.000 

IB 

0.000 

0 

■1 


0.000 

0.261 

0.125 

0.000 


-0.038 

1 



0.000 

0.781* 

0.375 

0.000 

i^ja 

-O.lll* 

3 



-0.125 





0.029 

0 

■■ 


-0.125 

^EallS 

0.125 

BUM 



1 


19 

-0.125 

0.1*09 


BiUM 


-0.085 

3 

■i 


-0.625 

-1.875 

0.000 




0 



-0.625 

-I.61I* 

0.125 

0.625 



1 

5 


-0.625 

-1.091 

0.375 

0.625 



3 




0.000 

0.000 



0.000 

0 




0.531 

0.125 



-0.035 

1 

0 



1.591* 

0.375 


-0.051 

-0.106 

3 



-0.125 

-0.525 

0.000 


0.008 

0.018 

0 

m 


-0.125 

0.006 

0.125 


-0.009 

-0.017 

1 


25 

-0.125 

1.069 

0.375 


-0.0U3 

-0.088 

3 

H 


-0.625 

1 

-2.625 

0.000 

IB 

0.01*2 

0.090 

0 



-0.625 

-2.091* 

0.125 


0.025 

0.055 

1 

5 


-0.625 

-1.031 

0.375 

IB 


-0.016 

3 




C-1* 











































































Table C-2. Coefficients of Errors dB, dH, dr in the Equations for the 
Errors d&z of Height Differences Az, and for Errors dAy of 
Distances Ay in a VOIR Stereo Arrangement of Figiire C-1. 


Error Coefficients 


Base 

dB 

tieight 

dH 

Range dr 


Ay 

Az 

Ay 

Az 

Ay 

Az 

Az 

Ay 

y 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0 



0,000 

0.021 

0.025 

0.000 

-0.001 

-0.001 

1 

0 


0.000 

0.063 

0.075 

0.000 

-0.003 

-0.003 

3 



-0.025 

-O.OU6 

0.000 

0.025 

0.001 

0.001 

0 

n 


-0.02,> 

-0.025 

0.025 

0.025 

0.000 

0.000 
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